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Universal features of phonon transport in nanowires with correlated surface 
roughness 
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The ultralow thermal conductivity ft observed experimentally in intentionally roughened silicon nanowires 
(SiNWs) is reproduced in phonon Monte Carlo simulations with exponentially correlated real-space rough 
surfaces similar to measurement [J. Lim, K. Hippalgaonkar, S. C. Andrews, A. Majumdar, and P. Yang, 
Nano Lett. 12, 2475 (2012)]. Universal features of thermal transport are revealed by presenting k as a 
function of the normalized geometric mean free path A (0 < A < 1); the diffusive (Casimir) limit corresponds 
to A = 1/2. ft vs A is exponential at low-to-mo derate roughness (high A), where internal scattering randomly 
interrupts phonon bouncing across the SiNW, and linear at high roughness (low A), where multiple scattering 
events at the same surface results in ultralow, amorphous-limit thermal conductivity. 


Heat flow through semiconductor nanostructures is 
governed by phonons, the quanta of lattice waves. 1 
Nanoscale phonon transport is presently a very active 
field of research, 2,3 straddling basic science inquiry 4 6 
and applications in electronics, 7-9 optoelectronics, 10 12 
and thermoelectrics. 13-15 Aided by considerable advances 
in measurement, 16,17 the design space for thermal trans¬ 
port has been greatly expanded by the use of low¬ 
dimensional and nanostructured semiconductors, 18,19 
which range from quasiballistic graphene-based sys¬ 
tems with superior heat conduction 6,20,21 to struc¬ 
tures in which rough surfaces or interfaces lead to 
low and anisotropic thermal conductivity, 22-24 such as 
superlattices, 22,23,25,26 nanomembranes, 27-31 , nanowires 
(NWs), 14,32-35 and structures with nanodots. 36,37 

However, theoretical understanding of phonon dynam¬ 
ics in nanostructures with a significant degree of dis¬ 
order is far from complete. A prominent open prob¬ 
lem is the unexpectedly low thermal conductivity, ft, 
of very rough silicon nanowires (SiNW). 14,34,35 Earlier 
measurements on vapor-liquid-solid (VLS)-grown NWs 
showed ft an order of magnitude lower than the bulk 
value. 32 These results could largely be explained within 
the diffusive-transport framework and a simple model 
of phonon interaction with surfaces: 1,38 the rough sur¬ 
face is described by a specularity parameter p, which 
is the probability that a phonon would reflect specu¬ 
larly upon impact. The diffusive or Casimir limit cor¬ 
responds to complete momentum randomization at the 
surface (p = 0). Indeed, thermal transport in many 
relatively smooth NWs is well described using the re¬ 
laxation time approximation and the model of partially 
or completely diffuse surface scattering. 3,32,39-41 In con¬ 
trast, the ft’s measured on electrolessly etched 14 or in¬ 
tentionally roughened VLS-grown 34 NWs are far be¬ 
low the Casimir limit and have values similar to amor¬ 
phous materials. Several groups have calculated the 
ft of SiNWs in the presence of periodic or rough sur¬ 
face features by molecular dynamics, 42-44 elastodynam- 
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FIG. 1. Example surface profiles with exponential (left) and 
Gaussian (right) autocorrelation functions. They have the 
same rms roughness (2 nm) and correlation length (25 nm, 
chosen large for visual appeal), and were created from the 
same set of random numbers. 


ics or hydrodynamics, 45-48 solving the phonon Boltz¬ 
mann transport equation (BTE), 33,40,41,49 53 or within 
coherent transport formalisms. 54 58 The unusually low 
thermal conductivity has been linked to diverse mech¬ 
anisms, such as the need to include full dispersion, 41 
different surface-scattering rates for phonons of different 
energies, 56,59 a surrounding native oxide that might ran¬ 
domize phonon phase, 3,60 internal defects, 35,44 and mul¬ 
tiple coherent backscattering events from highly corre¬ 
lated surfaces. 57,58 Recently, Lim et a/. 34 performed a 
systematic experimental study of the surface condition of 
intentionally roughened VLS-grown NWs, revealing cor¬ 
relation lengths smaller than previously assumed 57 and 
underscoring the importance of surface-roughness profile 
in nanoscale thermal transport. 

In this letter, we show that the ultralow measured ther¬ 
mal conductivities in intentionally roughened SiNWs can 
be reproduced in phonon Monte Carlo simulation with 
exponentially correlated real-space rough surfaces similar 
to experiment. 34 We introduce the normalized geomet¬ 
ric mean free path A, a dimensionless quantity propor¬ 
tional to the NW volume-to-surface ratio, which serves 
as a universal quantifier of surface roughness; A = 1/2 
corresponds to the Casimir limit, ft (A) is exponential 
at low roughness (high A), where internal scattering 
competes with roughness scattering, and linear at pro¬ 
nounced roughness (low A), where multiple successive 
scattering events from the same surface dominate and 
result in ultralow, amorphous-limit ft. 

We consider long SiNWs with square cross-sections 
and rough surfaces characterized a given autocor- 
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relation function (ACF). If h(r) is the height of 
surface S at point r, then the ACF is given by 
C (r) = ff s h(r') h (r' + r) cPf' . Most previous 

studies that included correlated roughness used ei¬ 
ther Gaussian 45 ’ 47,57,59 ’ 61 or Gaussian-based 46 ACFs. 
However, measurements of Si/Si02 interfaces found 
exponential 62 or exponential/Gaussian hybrid ACFs. 63 
Hybrid ACFs have also been found on SiNWs. 34,58 There¬ 
fore, we consider surfaces with both exponential [C e (r) = 
A 2 e - l r l/£] and Gaussian [C g (r) = A 2 e~^ 2 ^ 2 ] ACFs (ex¬ 
ponential and Gaussian surfaces, for brevity), where A 
denotes the rms roughness and £ is the correlation length. 
In the simulation, each NW surface has a random rough¬ 
ness profile generated numerically according to a given 
ACF, A, and £. 53,64 Figure 1 illustrates the differences 
between exponential and Gaussian surfaces that have the 
same A and £. The two surfaces have similar large-scale 
features, but the exponential surface has considerably 
more small-scale roughness. 

The SiNWs we consider have widths 20-70 nm and are 
approximately 2 fim long, much longer than the room- 
temperature bulk phonon mean free path, 65 so trans¬ 
port is scattering-limited (diffusive) and described by 
the phonon BTE. 1,66 Phonon Monte Carlo (PMC) is 
an efficient stochastic technique for the solution to the 
Bte, 28,53,67-69 which can incorporate real-space rough¬ 
ness and simulate the SiNWs of experimentally relevant 
sizes. In PMC, a large ensemble of numerical phonons 
is tracked over time as they fly freely and undergo scat¬ 
tering from the real-space surface roughness (upon hit¬ 
ting a point on the rough surface, the phonon reflects 
specularly) or from internal scattering mechanisms (nor¬ 
mal, umklapp, and isotope scattering; we used the rates 
from Morelli et al. 70 ). The simulation includes longitu¬ 
dinal and transverse acoustic phonons; optical phonons 
carry little heat owing to their low group velocity, but 
they affect the scattering rates of acoustic phonons and 
this influence is included in the above rates. Bulk dis¬ 
persions are reasonable for NWs thicker than a few tens 
of nanometers 39-41 and, for simplicity, we use isotropic 
relationships for each acoustic branch. 71 Planar black- 
body contacts, whose temperatures differ by 20 K, are 
at the two ends of the wire. 69 With the contact tem¬ 
peratures given, the simulation runs until a steady state 
is reached, as evidenced by a linear temperature pro¬ 
file (constant temperature gradient VT) and a constant 
power flux along the NW. k is obtained from Fourier’s 
law, <f> = —kX7T. Implementation details can be found 
in Refs. 53 and 72. 

Figure 2 presents the results of PMC simulations for 
a large ensemble of 70-nm-wide SiNWs, akin to those 
measured by Lim et al , 34 The rough surfaces were gener¬ 
ated according to exponential or Gaussian ACFs, with a 
broad range of roughness parameters, A=0-5.5 nm and 
£=3-16 nm; this range includes experimental values. 34 
As expected, hi decreases with increasing A and increases 
with increasing £ (the surface “looks” smoother with a 
higher £. 73 ) Exponential surfaces have more small-scale 
roughness, and consequently lower /Ps, because they are 
more effective at scattering short-wave-length phonons 



FIG. 2. Room-temperature thermal conductivity k versus rms 
roughness A and correlation length £ for 70-nm-wide SiNWs 
with exponential (left panel) and Gaussian (right panel) au¬ 
tocorrelation functions. The Casimir limit for SiNWs of this 
width, approximately 42.7 W/m-K, is shown with a dashed 
white line. 



A (nm) 

FIG. 3. Thermal conductivity k versus rms roughness A for 
70-nm-wide SiNWs with different correlation lengths £ and 
autocorrelation functions of the surface roughness. Values 
from the PMC calculation with real-space exponential (solid 
curves) and Gaussian (dashed curves) surfaces are presented 
for the correlation lengths of 5 nm (red), 9 nm (green), and 
13 nm (blue). Thermal conductivity calculated by PMC 
with a momentum-dependent specularity-parameter model 
for phonon surface-roughness scattering (black solid curve) 
tends to the Casimir limit. The symbols represent experi¬ 
mental values from the supplement to Lim et al ., 34 obtained 
on intentionally roughened SiNWs of diameters 67-84 nm, and 
with correlation lengths 8-10 nm (green circles) and 11-15 nm 
(blue triangles) extracted from exponential fits to experimen¬ 
tally obtained ACFs. 


than their Gaussian counterparts. Furthermore, hi cal¬ 
culated with real-space roughness is below the Casimir 
limit (about 42.7 W/m-K for 70-nm NWs and shown in 
a white dashed line in Fig. 2) over a wide range of A 
and £ for both exponential and Gaussian surfaces. The 
Casimir limit was obtained from PMC, wherein the NW 
surfaces are flat but each phonon’s momentum is com¬ 
pletely randomized whenever it hits a surface. 

Figure 3 further illustrates the importance of account¬ 
ing for real-space correlated roughness. The specularity 
parameter 1,73 p is a wave concept that can be derived 
based on basic diffraction theory in the limit of no corre¬ 
lation as p(q) = exp (—4A 2 ^), where q± is the compo¬ 
nent of the wave vector normal to the surface. 30,73 Fig¬ 
ure 3 shows hi versus A for the specularity-parameter 
model (£ = 0) and for exponential and Gaussian surfaces 
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with different values of £ (5, 9, and 13 nm). The ex¬ 
ponential curves and the specularity model agree up to 
about A ^ 0.5 nm, which is on the order of the aver¬ 
age phonon wave length, after which the specular curve 
saturates at the Casimir limit (p = 0). Moreover, our 
calculation results for exponential surfaces are very close 
to experiments on intentionally roughened NWs of sim¬ 
ilar diameters and surface roughness parameters. 34 The 
slight discrepancy might stem from isotropic rather than 
full dispersions, 56,72 the fact that the measured SiNWs 
appear to have a hybrid rather than purely exponential 
ACF, 34 and not accounting for the native oxide. 44 

Thermal conductivity includes the effects of both in¬ 
ternal and surface scattering. Focusing on surface scat¬ 
tering alone, we note that a NW can be considered an 
open cavity. In a closed cavity, the geometric mean free 
path (GMFP), A, is the average distance a phonon trav¬ 
els between successive surface scattering events in the 
absence of internal scattering and is computed simply as 
A = , where V is the cavity volume and S is its surface 

area. 74 Our wires have only small openings at the ends, 
so they resemble closed cavities; numerical computation 
of A in our NWs shows them to be almost identical to 
the closed-cavity values. A long square NW of length L 
has V = LW 2 (since the average height of rough sur¬ 
faces is zero) and S > AWL (the rougher the NW, the 
larger its surface area), thus A < W . Therefore, we de¬ 
fine the normalized geometric mean free path, A = X/W 
(0 < A < 1), which enables us to compare across different 
wire diameters. Smaller A means greater roughness. 

Figure 4 shows A as a function of A and £ for expo¬ 
nential and Gaussian correlation on 70-nm-wide SiNWs. 
Note that the A(A,£) surfaces look qualitatively similar 
to A,£) (Fig. 2); both k and A decrease as A increases 
and increase as £ increases. The bottom two panels in 
Fig. 4 reveal an important universality: the contour plots 
of A(A,£) for 35 and 70-nm-wide NWs are nearly identi¬ 
cal for a given ACF. Therefore, once we specify the type 
of correlation, the normalized geometric mean free path 
A can be considered as a near-universal quantifier of sur¬ 
face roughness scattering, encompassing W, A, and £. 
From this point on, we will concentrate on exponential 
surfaces, as they resemble experiment more closely. 34 

Figure 5(a) shows k as a function of A for SiNWs of 
width 20, 35, and 70 nm and with exponential surfaces. 
Each plot presents thermal conductivity data for sev¬ 
eral hundred different NWs with real-space roughness. 
There are a number of interesting features in this graph. 
First, as A 1 (the smooth-surface limit), all three 
curves converge to the bulk thermal conductivity, 75 as 
expected. This smooth-wire limit is well captured by 
treating each surface with an appropriate specularity pa¬ 
rameter, as seen in Fig. 3. Second, the Casimir limit 
values (17.3, 26.1, and 42.7 W/m • K for the 20, 35, and 
70-nm NWs, respectively), obtained from PMC simula¬ 
tions with smooth but momentum-randomizing surfaces, 
occurs at A ~ 1/2, i.e. A = W/2. Therefore, thermal con¬ 
ductivity below the Casimir limit corresponds to A < 1/2. 

In the inset to Fig. 5(a) we present thermal conduc¬ 
tivity normalized to width, k/W , versus A when inter- 
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FIG. 4. (Top row) A, the geometric mean free path normal¬ 
ized to the SiNW width, as a function of rms roughness A 
and correlation length £ for a 70-nm-wide SiNW whose rough 
boundary surfaces are characterized by (a) exponential and 
(b) Gaussian autocorrelation functions. (Bottom row) Con¬ 
tour plots of the normalized geometric mean free path A ver¬ 
sus A and £ for exponential (c) and Gaussian (d) boundary 
surfaces. Here, solid lines correspond to NWs of width 70 
nm, while dashed lines represent 35-nm-wide NWs. Consec¬ 
utively colored contours correspond to 0.05 increments in A. 
The color scale is the same as in the top row. 


nal scattering mechanisms have be turned off and only 
surface-roughness scattering remains. For A > 1/2, i.e. 
above the Casimir limit, the curves for the three different 
NW widths coincide: roughness-limited n/W is univer¬ 
sal as a function of A in this regime. With increasing 
roughening, however, once the Casimir limit has been 
surpassed (A < 1/2), the curves start to diverge from one 
another. Therefore, there is a fundamental difference in 
surface-limited transport above and below the Casimir 
limit, an important issue that we explore further below. 

The most prominent feature of Fig. 5(a) is that, for 
each NW thickness, there is a crossover in k(X). In the 
low-to-mo derate roughness regime, k( A) is exponential 
(In k ^ A, with a width-dependent slope), notably so 
for the thicker two NWs (35 and 70 nm). Considering 
that these are the results of semiclassical simulations, 
localization 76 stemming from coherent superposition of 
waves is out of the question. We believe the expla¬ 
nation is related to the phenomenon of variable-range 
hopping, 77,78 where there is an exponential dependence of 
conductivity of the characteristic hopping range. What 
happens here, for low-to-mo derate roughness, can be re¬ 
ferred to as variable-range bouncing : phonon bounces 
between opposite sides of the NW, while the range vari¬ 
ability stems from the path being interrupted by inter¬ 
nal scattering events, as depicted in the top-left-corner 
schematic in Fig. 5(a). 

Another interesting aspect of Fig. 5(a) is the k( A) be¬ 
havior of very rough NWs (small A). First off, the cross¬ 
over happens at smaller A for thicker NWs, because in 
thicker NWs the relative importance of internal scatter- 
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FIG. 5. (a) Room-temperature thermal conductivity ft as 

a function the normalized geometric mean free path A for 
SiNWs of width 70 nm (blue squares), 35 nm (green dia¬ 
monds), and 20 nm (red triangles), obtained from PMC simu¬ 
lation on a large ensemble of SiNWs with real-space exponen¬ 
tially correlated rough surfaces. (Inset) ft normalized to width 
W versus A without internal scattering for the same NWs as 
in the main panel. Vertical dashed line denotes A = 1/2. 
(b) Thermal conductivity ft as a function of the absolute geo¬ 
metric mean free path A for very rough SiNWs [low-A region 
from (a)]. The solid black line corresponds to a simple RTA 
solution with r = X/v± that yields a thermal conductivity of 
ft = AX, A « 2 x 10 9 W/m 2 K. The schematics in the top 
left corners depict phonon trajectories, interrupted by inter¬ 
nal scattering events, in NWs with (a) low-to-moderate and 
(b) pronounced roughness. 


ing at a given A is higher. The small-A region is also 
harder to reach in thicker NWs, like the 70-nm one, be¬ 
cause very high A would be required, outside of the range 
we focused on here. The ft (A) behavior in the small-A re¬ 
gion can be explained by multiple scattering events at the 
same boundary. In Fig. 5(b), we zoom in on the low-A re¬ 
gion from Fig. 5(a) and present ft as a function of the ab¬ 
solute GMFP, A. We note that the thermal conductivities 
reach very low values, of only a few W/m-K, such as those 
measured by Hochbaum et al. 1A We also see that the data 
for 20-nm and 35-nm NWs fall on top of each other, which 
is in keeping with the intuitive picture that multiple scat¬ 
tering events from the same boundary govern transport, 
so crossing the wire, and thus the wire width, becomes 
irrelevant. Indeed, the low-A region agrees very well with 
the simple relaxation-time approximation (RTA) expres¬ 


sion ft = (27r) 3 f d 3 qc(b, q)v^(b, q)r(b, q), where we 
put r{b,q) = X/v±(b, q). Here, c(b , q) is the heat capacity 
of mode b at wave number q , while v\\(b,q) and v±(b, q) 
are the velocity components along and across the wire, 
respectively. Upon simple integration, we obtain ft = AX, 
where A « 2 x 10 9 W/m 2 K at 300 K [Fig. 5(b)]. This A 
is single-surface dominated [see top-left-corner schematic 
in Fig. 5(b)] and bears essentially no dependence on the 
NW width. We note that Sadhu and Sinha 57 argued 
that coherent transport and multiple correlated scatter¬ 
ing events at a surface are key to low thermal conduc¬ 
tivities. Here, we obtain ultralow thermal conductivi¬ 
ties within a semiclassical (incoherent) transport model. 
While coherent transport is clearly not required for ul¬ 
tralow ft, multiple scattering processes from a surface are. 

In summary, we showed that thermal conductivity ft 
far below the Casimir limit can be readily obtained for 
NWs within the Boltzmann transport formalism, pro¬ 
vided that real-space rough surfaces with realistic ACFs 
are employed. We introduced the concept of the nor¬ 
malized geometric mean free path A, which encompasses 
NW width, roughness rms height, and correlation length. 
Thermal conductivities below the Casimir limit corre¬ 
spond to A < 1/2. ft(A) reveals universal signatures of 
the interplay between boundary-roughness and internal 
scattering: a) in the low-roughness, high-A region, ft (A) is 
exponential, owing to phonons bouncing across the NW 
and having trajectories randomly interrupted by internal 
scattering events; b) in the high-roughness, low-A region, 
multiple scattering events at a single interface govern 
transport, ft ^ A, and extremely low values of ft, close to 
the amorphous limit, are obtained. This work shows that 
pronounced roughness results in a fundamental, qualita¬ 
tive change to thermal transport in nanostructures. 

Finally, while the NWs with our numerically generated 
surfaces are technically not truly chaotic cavities, 79-82 
the real NWs likely are. 34 An individual phonon in a 
real SiNW may exhibit classically chaotic 83 or quantum- 
chaotic dynamical features, 84 depending on its wave 
length and whether coherence is preserved upon surface 
scattering. 3 A deeper understanding of thermal transport 
in very rough NWs means understanding how ensembles 
of phonons behave inside chaotic cavities. 85 
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